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Thermalization is investigated for the one-dimensional anisotropic antiferromagnetic Heisenberg 
model with dimerized nearest-neighbor interactions that break integrability. For this purpose the 
time evolution of local operator expectation values after an interacting quench is calculated directly 
with the Chebyshev polynomial expansion, and the deviation of the diagonal from the canonical 
thermal ensemble value is calculated for increasing system size for these operators. The spatial and 
spin symmetries of the Hamiltonian are taken into account to divide it into symmetry subsectors. 
The rate of thermalization is found to weaken with the dimerization parameter as the Hamiltonian 
evolves between two integrable limits, the non-dimerized and the fully dimerized where the chain 
breaks up into isolated dimers. This conclusion is supported by the distribution of the local oper¬ 
ator off-diagonal elements between the eigenstates of the Hamiltonian with respect to their energy 
difference, which determines the strength of temporal fluctuations. The off-diagonal elements have a 
low-energy peak for small dimerization which facilitates thermalization, and originates in the reduc¬ 
tion of spatial symmetry with respect to the non-dimerized limit. For increasing dimerization their 
distribution changes and develops a single low-energy maximum that relates to the fully dimerized 
limit and slows down thermalization. 

PACS numbers: 05.30.Ch,05.70.Ln,75.10.Pq,37.10.Jk 


I. INTRODUCTION 


In recent years the field of out-of-equilibrium phenom¬ 
ena in thermally isolated quantum many-body systems 
has attracted significant attention |ll-H4l|. A part of an 
isolated system feels the rest of the system acting as a 
bath, resulting in probability flowing in and out of the 
subsystem. This can lead to weakly fluctuating opera¬ 
tor expectation yalues in the subsystem in the long-time 
limit, which can be described by a statistical ensemble 
with a small number of parameters. The nature of the 
statistical ensemble depends on the integrability proper¬ 
ties of the Hamiltonian. Unlike the classical case, inte- 

a bility is not precisely defined in quantum mechanics 
. The prevailing idea is that in a quantum-integrable 
system the number of local integrals of motion is propor¬ 
tional to the constituents of the system, something not 
required at the classical level. A non-integrable Hamilto¬ 
nian is expected to thermalize according to the canonical 
thermal ensemble. On the other hand, the conserved 
quantities of an integrable Hamiltonian have to be taken 
into account in the generalized Gibbs ensemble that de¬ 
scribes correlations in the long-time limit @,01. It also 
has been recently shown that the quench action formal¬ 
ism correctly captures the long-time limit steady state 
Hamiltonians where the proximity to integrabil¬ 
ity is controlled by a parameter are of particular inter¬ 
est, since it is possible to interpolate between the inte¬ 
grable and non-integrable cases in a well-controlled fash¬ 
ion 


In the case of weak integrability breaking it has been 
shown that thermalization can be particularly slow. For 
the fermionic 23|and the spinless fermionic dimerized 
Hubbard model |24| . as well as the quantum Ising chain 


local operator expectation values first relax to a 


non-thermal quasi-stationary value in a prethermalized 
regime, before eventually relaxing to the canonical ther¬ 
mal ensemble value. Different scenarions of prerelaxation 
close to integrable points have been considered in [2^ . In 
the case of the anisotropic antiferromagnetic Heisenberg 
model (AAHM) with nearest and next-nearest neighbor 
interactions that break integrability, thermalization is 
also slow close to the integrable point [2lj. This was 
attributed to the appearance of a low-energy peak for 
local operator off-diagonal elements that facilitates ther¬ 
malization once integrability is broken, but which is only 
weak for small integrability breaking. For stronger inte¬ 
grability breaking the low-energy peak increases in mag¬ 
nitude and thermalization becomes faster. Altogether 
there are three different thermalization regimes when the 
next-nearest neighbor interaction is relatively moderate. 

In this paper we investigate thermalization for the 
AAHM for a spin chain with dimerized interactions, 
where the proximity to integrability is controlled by the 
dimerization strength. This is defined as the deviation 
of the nearest-neighbor exchange constants from the uni¬ 
form interactions limit. In the absence of dimerization 
the AAHM is integrable with the Bethe ansatz , how¬ 
ever an infinitesimal dimerization breaks integrability. At 
the opposite end of full dimerization the model is trivially 
integrable, with the chain breaking up in non-interacting 
dimers. Tuning the dimerization strength allows con¬ 
trol of the proximity to integrability, in analogy with the 
next-nearest neighbor interaction in the one-dimensional 
chain [2l| . To investigate the effect of integrability break¬ 
ing on the rate of thermalization we calculate the time 
evolution of local correlations after an interacting quench, 
and more specifically their deviation from their thermal 
ensemble value as a function of time. We also calculate 
the difference between the diagonal and thermal ensem¬ 
ble values for local operators for increasing system size. 
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This difference should go to zero in the thermodynamic 
limit (TDL) in the case of thermalization. In both calcu¬ 
lations the dependence on the dimerization parameter is 
examined. As in the case of the AAHM with next-nearest 
neighbor interactions, we invoke the distribution of the 
operator off-diagonal elements with respect to the energy 
difference to determine the relative strength of temporal 
fluctuations for varying dimerization. 

An infinitesimal dimerization breaks the integrability 
of the one-dimensional uniform AAHM and reduces its 
spatial symmetry by a factor of two by doubling the unit 
cell. This is qualitatively different with respect to the 
breaking of integrability by the addition of next-nearest 
neighbor interactions, which in contrast preserve the spa¬ 
tial symmetry of the Hamiltonian. In the latter case the 
operator off-diagonal elements that control the strength 
of temporal fluctuations of the operator’s expectation 
value change smoothly as integrability is broken, and the 
time evolution takes place in the same symmetry sub¬ 
sector as in the integrable case. On the other hand, 
dimerization mixes symmetry subsectors of the uniform 
AAHM due to the spatial symmetry reduction. This 
means that even an infinitesimal dimerization changes 
qualitatively the form of the operator off-diagonal ele¬ 
ments, as levels very close in energy immediately develop 
considerable overlaps. Such strong off-diagonal matrix 
elements between states close in energy are absent in the 
non-dimerized limit, and are responsible for the fastest 
thermalization occuring for very weak dimerization. In 
contrast, the AAHM with next-nearest neighbor inter¬ 
actions requires a stronger finite integrability breaking 
to generate significant overlaps between eigenstates not 
differing much in energy, and therefore to cause faster 
thermalization [2lj. 

In the opposite limit of full dimerization the AAHM 
is trivially integrable, with each non-zero bond connect¬ 
ing two spins to form an isolated dimer. For strong but 
not full dimerization the isolated dimers start to inter¬ 
act and form energy bands. Their interaction is however 
considerably weaker than intra-dimer exchange, conse¬ 
quently the probability flow between neighboring dimers 
is weak and thermalization is slow. This can also be seen 
in the off-diagonal elements of local operators between 
different eigenstates of the Hamiltonian. These opera¬ 
tors act on isolated or neighboring dimers of the fully 
dimerized limit, which means that they will have non¬ 
zero off-diagonal elements between eigenstates whose en¬ 
ergies differ at the level of one or two dimers and are 
therefore quite close. 

According to the above, the breaking of integrability of 
the dimerized AAHM is qualitatively different in the two 
extremes of weak and strong dimerization, as in the for¬ 
mer case it assists thermalization in a much more efficient 
way due to the immediate reduction of spatial symme¬ 
try. In the latter case the evolution away from the fully 
dimerized limit is smooth and efficient thermalization set¬ 
tles in more slowly, as in the case of the AAHM with 
next-nearest neighbor interactions. It is found that the 


time-dependent local operator expectation values evolve 
smoothly between the non-dimerized and the fully dimer¬ 
ized limit, meaning that thermalization is getting slower 
as the dimerization is getting stronger. Further evidence 
for this is provided from the operator off-diagonal ele¬ 
ments, which also evolve with the dimerization smoothly 
between the two extreme limits, and also from the dif¬ 
ference between the diagonal and thermal ensemble val¬ 
ues. The result is that for moderate dimerization ther¬ 
malization is already quite slower. Here the interaction 
term of the AAHM is chosen comparable in magnitude 
to the tunneling term, with the interaction also stronger 
than the dimerization. Similar results for considerable 
dimerization have been found when the interaction term 
is weaker than the tunneling term [^ . where the slow 
approach to thermalization was associated with the emer¬ 
gence of the prethermalized regime. 

The plan of this paper is as follows: Section m intro¬ 
duces the model and the methods to be used. Section Hill 
discusses the time evolution of an operator expectation 
value after a quench and the requirements for thermal¬ 
ization. Section m includes the results, namely time- 
dependent local operator expectation values and the scal¬ 
ing of the difference between the diagonal and canonical 
thermal ensemble values for an interacting quench. It 
also presents the distribution of the strength of local op¬ 
erator off-diagonal terms with energy difference. Finally 
Sec. rvl presents the conclusions. 


II. MODEL AND METHOD 

The Hamiltonian of the dimerized AAHM model for a 
chain is 

N 

H = J^[(I - -f Asfsf+JI) 

i=l 

with N the number of spins. The spin magnitude 
s = 1/2 and the boundary conditions are periodic so that 
sn+i = Si . In contrast to the non-dimerized case the unit 
cell contains two spins. In order for the ground state be¬ 
fore the quench and consequently the time-evolved wave- 
function to be included in the one-dimensional even ir¬ 
reducible representation [2^ , the strong and weak bonds 
must both appear an even number of times, therefore 
N is taken to be a multiple of four. The dimerization 
parameter 5 is taken to be positive and breaks the in¬ 
tegrability of Hamiltonian o once it is non-zero, with 
the exchange constants alternating as J(l±(5). Hamilto¬ 
nian ([T]) is trivially integrable in the fully dimerized limit 
(5 = 1. J is taken to be 1 from now on, defining the unit 
of energy. 

To calculate diagonal and thermal ensemble values 
Hamiltonian m has to be fully diagonalized, and in or¬ 
der to do so its symmetries are exploited [29l - [^ . The 
full symmetry group of the Hamiltonian is the product of 
the spatial and spin symmetry groups. The spatial group 
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is the dihedral Dn_ group [l^, while in spin space the 
Hamiltonian is symmetric under inversion when = 0. 

To calculate directly the time evolution of opera¬ 
tors the Chebyshev polynomial expansion is employed 
[allai. It is more efficient in comparison with exact 
diagonalization for this purpose, making it possible to 
consider larger system sizes. In addition, the block- 
diagonalization of Hamiltonian ([1]) according to the ir¬ 
reducible representations of its total symmetry group al¬ 
lows to time-evolve the initial wavefunction separately 
within each representation, increasing the maximum sys¬ 
tem size that can be considered even more. If the initial 
wavefunction is confined to a single symmetry subsector, 
it evolves in time only in this subsector and furthermore 
the diagonal ensemble has non-zero contributions only 
from expectation values of this subsector. This is the 
case for the interacting quench considered in this paper. 
The initial state is the ground state of Hamiltonian (ED 
before the quench. For the initial state different values of 
Ao are chosen for different 5 so that the post-quench ef¬ 
fective temperature is equal to P = 1.3500 in every case. 
After the quench A = 1.1 is taken. 

In the absence of dimerization {6 = 0) Hamiltonian ([T]) 
is integrable via the Bethe ansatz [23| , and there is quasi- 
long range order in the ground state. Once S ^ 0 a gap 
opens and the ground state has Neel order [^1^. Here 
we consider relatively small <5 < 0.3, and investigate how 
integrability breaking affects thermalization by consider¬ 
ing the nearest-neighbor operators + and s|s|, 

on the strong (j = 2i + 1, k = 2i + 2) and weak (j = 2i, 
k = 2i + 1) bonds. We pick A = 1.1 for the anisotropy 
parameter after the quench to rule out mixing of differ¬ 
ent S sectors, which can occur in the spin isotropic SU(2) 
case when one works in the basis. 


III. TIME EVOLUTION 

The initial wavefunction |d>(t = 0)) = |d>(0)) is the 
ground state of the Hamiltonian before the quench. If the 
post-quench eigenvalues and eigenvectors of Hamiltonian 
ED are En and |it is for the wavefunction |d>(t)) 
after the quench 

|^(t)) =e-*^‘|^(0)) (2) 

n 

The coefficients C„ = (d>„|'I'(0)) do not depend on time. 
The time evolution takes place in the one-dimensional 
even irreducible representation (Sec. HID: therefore there 
are no degeneracies and an operator O at time t has the 
expectation value 

<d{t) > = + 

n 

Y, (3) 

m^n 


In the time average of Eq. m only the first term has 
non-zero contribution. The average is given by O = 
J2n |C'nP('I'„^|'I'n), defining the value of O in the diago¬ 
nal ensemble [7| . The wavefunction evolves in time in the 
non-degenerate subsector of the ground state before the 
quench, therefore only eigenstates of this subsector have 
non-zero coefficients. Thermalization requires the diago¬ 
nal and canonical thermal ensemble values to be equal in 
the TDL. According to Eq. ([3D the strength of fluctua¬ 
tions around the diagonal ensemble value is controlled by 
the off-diagonal terms ('I'„|0|'I'm) = Onm, m ^ n. Ther¬ 
malization also requires that fluctuations with respect to 
O are small in the TDL. To determine the canonical ther¬ 
mal ensemble pth describing the equilibrated long-time 
limit, the energy and the thermal average are set equal 
to each other 

{^mH\^{Q))=Tr{Hpth) (4) 

The numerical solution of this equation gives the tem¬ 
perature. 


IV. NEAREST-NEIGHBOR CORRELATIONS 

To investigate thermalization the time evolution of lo¬ 
cal operator expectation values is of interest, therefore 
nearest-neighbor interactions are considered, which re¬ 
late to a subsystem of two adjacent spins. The full cor¬ 
relation Si ■ Si+i is not calculated, as it commutes with 
Hamiltonian ED at the integrable points and for isotropic 
spin interactions A = 1. Instead correlations in the spin 
xy plane and along the spin z axis are considered. Fig. 
[T] shows the time evolution of the expectation value of 
the planar spin correlation S 2 i_|_iS 2 i _|_2 + S2i-i-i'52i-i-2 
tween spins interacting via the strong exchange coupling 
1-1-(5. More accurately, it shows the deviation of the time- 
dependent expectation value from the canonical thermal 
ensemble value. Since the chains considered are of fi¬ 
nite size, revivals of expectation values are eventually 
generated in later times. The results presented are for 
a V = 32 chain, which generates the infinite chain re¬ 
sults without any revivals for times at least up to j 
[m, [ 2 ^. This has also been vindicated by comparing 
the time-dependent expectation values as a function of 
V, as they agree to progressively longer times with in¬ 
creasing N. The deviation from zero in Fig. |T] increases 
with the dimerization, indicating that thermalization be¬ 
comes slower with 5. Further evidence of that comes from 
the fluctuations for higher dimerizations, which are still 
of considerable strength for the highest times. 

It is instructive to compare the case of the dimerized 
Hamiltonian ED with the integrability breaking effected 
by adding uniform couplings between next-nearest neigh¬ 
bors. In this case thermalization becomes faster with the 
strength of the integrability breaking term [2l|. This 
was explained by considering the energy distribution of 
the off-diagonal elements of the local operators, which 
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according to Eq. © control the strength of the fluc¬ 
tuations around the long-time limit. It was found that 
for very small next-nearest neighbor exchange the off- 
diagonal elements peak away from the origin. As the 
next-nearest neighbor exchange increases, a peak close 
to zero energy develops that becomes important in the 
long-time limit and facilitates thermalization. This is 
because the high-frequency fluctuations related to the 
higher-energy peak will eventually cancel each other out, 
and the low-energy peak will determine the fluctuations 
in the long-time limit, which will be slowly varying. Fig. 
[2]shows the same calculation for the off-diagonal elements 
of ^ 2 i+i^ 2 i +2 + ^ 2 i+i^ 2 i +2 post-quench Hamilto¬ 

nian. In the integrable limit 5 = 0 the Hamiltonian corre¬ 
sponds to the AAHM with zero next-nearest neighbor in¬ 
teractions, where the local operator off-diagonal elements 
are very small close to zero and have a maximum away 
from it For very small 5 = 0.001 their distribution 
is quite different. Matrix elements between states very 
close in energy have now values comparable but smaller 
in magnitude than the higher-energy peak. Even though 
there is no well separated low-energy peak as in the case 
of the next-nearest neighbor AAHM, the present distri¬ 
bution also points to thermalization in the long-time limit 
in accordance with Fig. [U where the low-frequency fluc¬ 
tuations will eventually overtake the high-frequency ones. 

The introduction of dimerization breaks the integra- 
bility of Hamiltonian (IT|) and alters significantly the dis¬ 
tribution of the off-diagonal elements. The time-evolved 
wavefunction belongs to a specific irreducible represen¬ 
tation, which is determined by the ground state before 
the quench. In the integrable limit 5 = 0, eigenstates 
close in energy that belong to the irreducible representa¬ 
tion that solely includes the time-evolved wavefunction 
have local operator off-diagonal elements close to zero 
(see discussion in Ref. [Hi). An infinitesimal 5 reduces 
the spatial symmetry of the Hamiltonian by a factor of 
two by doubling the unit cell, and the spatial symmetry 
group changes from Df^ to Dn [2^. The total number of 
irreducible representations is less for the group, and 
their eigenstates are combinations of eigenstates of differ¬ 
ent irreducible representations of the Dn group. There 
is no restriction on the strength of the off-diagonal ele¬ 
ments of local operators between eigenstates with similar 
energies that belong to different irreducible representa¬ 
tions of the Dn group. The appearance of significant 
local operator off-diagonal elements between eigenstates 
close in energy for infinitesimal dimerization is thus not 
a perturbative effect, but originates in the reduction of 
the Hamiltonian spatial symmetry, which also results in 
the breaking of integrability. This is seen in Fig. [5] for 
5 = 0.001 and 0.01. 


In the opposite limit of strong dimerization Hamilto¬ 
nian © is trivially integrable for 5 = 1, with each non¬ 
zero bond connecting two spins to form an isolated dimer. 
The eigenstates of the Hamiltonian are products of dimer 
eigenstates, with their eigenenergies being sums of dimer 
energies, which are — | and scaled with the intra-dimer 


exchange interaction. For strong dimerization 5 < 1 the 
isolated dimers start to interact and form narrow energy 
bands around the discrete levels of the fully dimerized 
limit. The inter-dimer interaction is however consider¬ 
ably weaker than the intra-dimer exchange, consequently 
the probability flow between neighboring dimers is weak, 
which points to slow thermalization. This is in agreement 
with Fig. [U where the deviation from zero is considerably 
stronger for larger 5 and also fluctuations are strong. In 
addition, local operators act on approximately isolated 
dimers or weakly connected neighboring dimers, there¬ 
fore they will only have non-zero off-diagonal elements 
between eigenstates whose energies differ at the level of 
one or two dimers, and are therefore close and belong 
to the same energy band. This is seen in Fig. [2J where 
for the higher 5’s the off-diagonal elements have a maxi¬ 
mum for an energy difference close to zero, and decrease 
with increasing energy difference. This off-diagonal ele¬ 
ments distribution is characteristic of a strongly dimer¬ 
ized Hamiltonian, with the strong low-frequency fluctu¬ 
ations leading to slow thermalization. Fig. [5] also shows 
that the evolution of the off-diagonal elements distribu¬ 
tion of -I- S 2 i-i-i® 2 i-i -2 eigenstate basis is 

smooth as a function of 5 between the weak and strong 
dimerization limits. This points to slower thermalization 
for increasing dimerization, in agreement with Fig. [T] 

Further evidence for the dependence of the thermal¬ 
ization rate on 5 comes from the scaling of the difference 
between the diagonal and the canonical thermal ensemble 
values with N. This difference scaled with the canoni¬ 
cal thermal ensemble value is shown in Fig. [3Ka). For 
the two smallest 5 values the dependence on ^ appears 
very close to linear. If a straight line is fitted through the 
data, the intercept is equal to 3.6 x 10“^ for 5 = 0.001 and 
—7.2 X 10“^ for 5 = 0.01. This shows that chains of max¬ 
imum size A = 24 point to thermalization in the TDL 
for weak dimerization, with thermalization being faster 
for smaller 5. For the next highest 5 = 0.05 the line is 
significantly curved, and even a fit through the last two 
points does not pass close from the origin. This points 
to slower thermalization in agreement with Fig. [2 where 
for 5 = 0.05 the low and the high-energy peak become 
approximately of the same strength. For even higher 5 
the lines are also curved, and also the distribution of the 
off-diagonal elements has its highest weight close to zero, 
both indicating even slower equilibration to a long-time 
limit as the Hamiltonian gets even closer to the strongly 
dimerized limit. Due to the lack of quicker thermaliza¬ 
tion one would have to go to higher N to deduce the 
dependence of the ensemble difference on A. 

Fig. ISKb) shows the corresponding result for correla¬ 
tion function S 2 i+is|j_|_ 2 . A similar conclusion can be 
drawn. The intercepts for a straight line fit through 
the data for 5 = 0.001 and 0.1 are 0.011 and 0.016 re¬ 
spectively, indicating faster thermalization for smaller 5, 
while for higher 5 the dependence on A becomes again 
non-linear, at least for the available sizes. In Fig. 0] 
the difference of the time-dependent expectation value of 
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^ 2 i+i^ 2 i +2 from its thermal ensemble value is plotted as a 
function of time. <5 = 0.001 is not the closest curve to zero 
for the longest available time, however for higher S the de¬ 
viation has already changed sign and this is also expected 
for the deviations for 6 closer to 0.001. This conclusion 
is also supported by the scaling data of Fig. El^b), which 
resemble a straight line best for weak 6, and also from the 
distribution of the off-diagonal elements, which is shown 
in Fig. [5] and is similar to Fig. Again the results 
for correlation function S 2 j_|_]^S 2 j _|_2 support the thermal- 
ization scenario where stronger dimerization slows down 
thermalization. 

The fluctuations in Fig. |4]are the strongest for higher 
S. Fluctuations for larger <5 are expected to be significant, 
judging also from the distributions of Fig. [5l The strong 
fluctuations of the difference around zero for the higher S 
lead to a relative difference in Fig. |3ljb) that approaches 
zero fast with system size, even though the subsystem 
has not thermalized yet. This is something that does 
not occur for S 2 *-Hi'S 2 i-h 2 + sli+A +2 or in the AA^ 
with nearest and next-nearest neighbor interactions [21| . 
The reason is that in these cases the difference does not 
change sign for later times. 

Similar conclusions for slower thermalization with 
increasing 5 are drawn from the correlation function 
^ 2 i^ 2 i+i + ^ 2 i^ 2 i+i between nearest-neighbor spins inter¬ 
acting via the weak exchange coupling 1 — (5. Fig. El^a) 
shows the scaling data, where for the smallest S values 
0.001 and 0.01 a straight line fit passes very close from the 
origin. For higher S this is not the case, and the depen¬ 
dence on A deviates from a straight line, at least for the 
available N. Fig. [7] shows the time evolution of the ex¬ 
pectation value, which is similar to the one of S 2 i-i- 1 * 21 - 1-2 
and also points to slower thermalization with increasing 
dimerization. In Fig. [3 it is expected that for small <5 
the deviations will also change sign for longer times, as 
is also the expectation for Fig. ID The same thermal¬ 
ization pattern is suggested by correlation S2^S2i+l from 
the results shown in Figs. [SJb) and H) For both corre¬ 
lation functions on the weak bonds the thermalization 


predictions are supported by the distributions of the off- 
diagonal elements of the operators, which are not shown 
but are very similar to the ones of Figs. [2]and[5l 


V. CONCLUSIONS 

It is generally expected that a non-integrable Hamil¬ 
tonian thermalizes according to the canonical thermal 
ensemble. In this paper we considered a Hamiltonian 
that interpolates between two integrable points by vary¬ 
ing the degree of dimerization of its nearest-neighbor in¬ 
teractions. It was found that thermalization slows down 
with the dimerization strength. This was shown by three 
different types of calculations. First the direct time evo¬ 
lution of expectation values of local operators with time 
was computed, as well as the difference between the di¬ 
agonal and thermal ensemble values with increasing sys¬ 
tem size for these operators. Finally the distribution of 
local operator off-diagonal elements as a function of en¬ 
ergy, which control the temporal fluctuations, was calcu¬ 
lated. These calculations show that the reduction of spa¬ 
tial symmetry leads to the quickest thermalization for in¬ 
finitesimal dimerization. This contrasts the AAHM with 
nearest and integrability breaking next-nearest neighbor 
interactions, where a next-nearest neighbor interaction 
of si gnif icant strength is required to speed up thermaliza¬ 
tion [2l| . In the fully dimerized limit the chain breaks up 
into isolated dimers which do not exchange probability, 
therefore there is no thermalization and the off-diagonal 
elements peak close to zero. The calculated quantities 
evolve smoothly between the non-dimerized and the fully 
dimerized limit and show how increasing dimerization 
slows down thermalization. A dimerized Hamiltonian is 
also important as it can be considered the first step on 
the way to a random Hamiltonian, where randomness can 
lead to a complete lack of thermalization [s^ . 

The author is very thankful to A. Lazarides. 
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FIG. 1: (Color online) Time evolution of the expectation value 
of O = S 2 i+iS 2 i +2 + ^ 2 i+i^ 2 i +2 with respect to its thermal 
value < O >th for different values of 5 with A = 1.1. Spins 
S 2 i+i and S 2 i +2 interact via a strong bond 1 + 5. For the initial 
state different values of Aq are chosen for different S so that 
the post-quench effective temperature is equal to /? = 1.3500 
for every S. Starting with the weakest and ending with the 
strongest S the values of Aq are 88.94, 88.236, 82.95, 74.152, 
64.1131, 50.124, 33.4604, 23.7088, 18.0157, and 14.5974. 



FIG. 2: (Color online) Distribution of matrix elements 
{s 2 i+iS 2 i +2 + S 2 i+i^ 2 i+ 2 )r>.m as a function of the energy dif¬ 
ference AEnm = En — Em for different values of S with 
A = 1.1. What is plotted is the coarse grained average of 
l(•S2^+l■52^+2 + S 2 i+i^ 2 i+ 2 )nm \ divided with its maximum value 
as a function of the reduced average energy difference for 100 
bins with an equal number of points. Spins S 2 i+i and S 2 i +2 
interact via a strong bond 1 + 5. 
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FIG. 3: (Color online) Scaling of the difference AOrei = 
between the diagonal and canonical thermal ensem- 

\^}th 

ble values divided by the canonical thermal ensemble value 
as a function of inverse length i for correlation function (a) 

o = S2i+iS2i+2 + S 2 i+iS 2 i +2 ^nd (b) 6 = sii+isii +2 for differ¬ 
ent values of 5 with A = 1.1. Spins S 2 i+i and S 2 i +2 interact 
via a strong bond 1-1-5. The maximum chain has N = 24. 
The post-quench effective temperature is equal to /3 = 1.3500 
for every 5, and the values of Aq are given in the caption of 

Fig.m 
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FIG. 4: (Color online) Time evolution of the expectation value 
of O = S 2 i+iS 2 i +2 with respect to its thermal value < O >th 
for different values of S with A = 1.1. Spins S 2 i+i and S 2 i +2 
interact via a strong bond 1-1-5. The post-quench effective 
temperature is equal to /3 = 1.3500 for every 5, and the values 
of Aq are given in the caption of Fig. [T] 
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FIG. 5; (Color online) Distribution of matrix elements 
is 2 i+iS 2 i+ 2 )nm US a function of the energy difference AEnm = 
En — Em for different values of 5 with A = 1.1. What is plot¬ 
ted is the coarse grained average of \(s 2 i+iS 2 i+ 2 )nm\ divided 
with its maximum value as a function of the reduced average 
energy difference for 100 bins with an equal number of points. 
Spins S2i-i-i and S2i+2 interact via a strong bond 1-1-5. 



FIG. 6: (Color online) Scaling of the difference AOrei = 
between the diagonal and canonical thermal ensem- 

\^}th 

ble values divided by the canonical thermal ensemble value 
as a function of inverse length i for correlation functions (a) 

O = sfisfi+i + S 2 iS 2 i-)-i and (b) O = S 2 iS 2 i+i for different 
values of 5 with A = 1.1. Spins S 2 i and S 2 i+i interact via 
a weak bond 1 — 5. The maximum chain has N — 24. The 
post-quench effective temperature is equal to /? = 1.3500 for 
every 5, and the values of Aq are given in the caption of Fig. 

ra 
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FIG. 7: (Color online) Time evolution of the expectation value 
of O = S 2 iS 2 i+i + S 2 iS 2 i+i with respect to its thermal value 
< O >th for different values of 5 with A = 1.1. Spins S 2 i 
and S 2 i+i interact via a weak bond 1 — <5. The labels for the 
different lines are the same as in Figs. [T] and [4] The post- 
quench effective temperature is equal to /3 = 1.3500 for every 
S, and the values of Ao are given in the caption of Fig. [T] 



FIG. 8: (Color online) Time evolution of the expectation value 
of O = S 2 iS 2 i+i with respect to its thermal value < O >th 
for different values of 5 with A = 1.1. Spins S2i and S2i+i 
interact via a weak bond 1 — 5. The labels for the different 
lines are the same as in Figs. [T] and [H The post-quench 
effective temperature is equal to /3 = 1.3500 for every 5, and 
the values of Ao are given in the caption of Fig. [1] 






















